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A dispersion equation for the interface waves is derived for the interface of two cubic crystals in 
£SJ \ the plane perpendicular to [001]. A reasonable hypothesis is made about the total number of the 

acoustic modes. Due to this hypothesis the number is 64, but not all of the modes have physical 
meaning of the interface waves. The rules have been worked out to select physical branches among 
all 64 roots of dispersion equation. The physical meaning of leaky interface waves is discussed. 
The calculations have been made for the interface Alo.3Gao.7As/GaAs. In this case all physical 
interface modes have been shown to be leaky. The velocities of the interface waves are calculated 
C^l . as a function of an angle in the plane of interface. The results support a recent interpretation of a 

new type oscillations of magnetoresistance as a resonant scattering of two-dimensional electron gas 
by the leaky interface phonons. 
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I. INTRODUCTION 



The propagation of acoustic surface and interface waves has attracted a significant attention over the last decades. 
The concept of the surface waves goes back to the famous paper by Lord RayleigbJj. The interface wave is a simple 
I generalization of the surface wave, when the second medium is not vacuum, and the wave propagates along the 
boundary between two media. The theoretical study of the interface waves was initiated by Stoneleyn who considered 
O ' the case of two isotropic solids. \—\n 

In anisotropic materials the interface waves between hexagonal crystalsEHo have been theoretically studied in 
sufficient details. Relatively little known about effects of crystalline anisotropy when the interface is formed by 
cubic crystals. To the best of our knowledge the only nMSBepcal search for true interface wave velocities for several 
J> ■ combinations of the materials has been performed so farc^ES. 

Both surfaca-and interface waves were initially studied in the context of seismological waves propagating in thc 
Earth's crustlTQ. Later on these waves have been studied experimentally in semiconductors by the light scattering!™. 

The earlier theoretical studies by Lord Rayleigh and Stoneley ( see alsoa) prescribe to consider only those roots of 
the secular equations which give an exponential decay of the surface wave in the medium under the surface and an 
exponential decay of the interface wave in both media away from the interface. 

Probably, Phinneyfl was the first to consider the so called "leaky" or "pseudo" waves which do not obey this 
prescription. Surface leaky waves have been widely studied theoretically for both isotropic and anisotropic crystalline 
materials (see brilliant review by Maradudk£3 and references therepj). To the best of our knowledge leaky interface 
waves were studied only for isotropic caseOEl and hexagonal crystalail 

The interest to the leaky interface waves is stimulated by the fact that true interface waves exist inside a very 
narrow range of the parameters. Therefore in general case interface waves are leaky. This is not the case for the 
surface waves where the true non-leaky mode always exists in a wide range of parametersEj. However, the dispersion 
£j , equation for surface waves also has several roots which give leaky solutions. 

The structures with two-dimensional electron gas (2DEG), like heterostructures or quantum wells provide another 
source of an interest, to the interface waves. The study of the interaction of 2DEG with surface waves has been 
investigated long agcOEj. If 2DEG is far from the surface, the electrons may interact with interface waves. Say, the 
electrons may be scattered by thermally excited interface waves. This scattering should not be less than scattering 
by bulk phonons, since in the vicinity of the-Lnterface the three-dimensional densities of the bulk and the interface 
phonons are of the same order. In the papeiuJ we explained the novel oscillations of magnetoresistance, observed in 
high- mobility 2DEG in GaAs-AlGaAs heterostructures, by magneto-phonon resonance originating from interaction of 
the 2DEG with thermally excited leaky interface acoustic phonon modes. 

The primarily goal of this paper is calculation of the interface waves for Alo.3Gao.7As/GaAs interface on the basal 
(001) face. This is exactly the interface used in the paperEj. We have shown that all interface waves in this case are 
leaky. 

To this end we have derived analytically the secular equation for phase velocity v of the waves at the interface 
between two cubic crystals. We have discussed the selection rules for the modes and have given a novel general 
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qualitative picture of the leaky interface waves. In this picture we consider the conservation of energy and show 
that the amplitude of the wave never becomes infinite if the problem is properly formulated. We show that at some 
conditions leaky waves does not differ substantially from the true waves. Finally we have obtained the numerical 
results, which were partially used in Ref.Ej. 

The paper is organized as follows. The basis of the method is outlined in section [n]. In the third section we discuss 
general properties of the secular equation, the selection rules for its solutions, and the physical meaning of leaky 
waves. The numerical results and discussion are presented in section IV. Finally, some auxiliary technical material 
regarding calculations is given in the Appendices. 



II. GENERAL FORMULATION 



Within the framework of the linear theory of elasticity the equations of motion of the infinite medium are 

du, dan 



P 



dt 2 dxj 



1,2,3, 



(1) 



where p is the mass density of the medium, m(r,t) is Cartesian component of the displacement of the media at the 
point r at time t, and aij(r,t) is the stress tensor. The latter is given by Hooke's law 



dm 



Gil ^ijkl 5 

dxi 



where Xijki is symmetrical forth rank tensor. In cubic crystal stress tensor can be conveniently written as 



crij = Ci 2 (div u)5ij + C 4 



dui duj 
dxj dxi 



n dUl x 
D dx~ 5 ^ 



(2) 



(3) 



where summation over i and j is not assumed in the last term. Here D is the anisotropy parameter: 

D = C n - Cx2 - 2C 44 . (4) 

For isotropic medium D — 0. 

In the following analysis we consider a system formed by two semi- infinite cubic crystals. Elastic constants and 
density related to the lower part will be denoted by prime symbol (see Fig. 1). The interface is supposed to be on 
(001) cut and perpendicular to x 3 axis. 
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FIG. 1. Structure for the study of interface acoustic waves. 
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The equations of motion ([!]) have to be supplemented by the boundary conditions on the interface, expressing 
continuity of the displacement and normal components of the stress tensor: 

Ui = u\ \ X3=0 , i = 1, 2, 3, (5) 
Oi3 = 0*3 1x3=0, i = 1,2,3. (6) 

The homogeneous plane waves (bulk phonons) are the simplest solutions of the wave equation in one infinite medium. 
They are 

u^(r, t) — exp (ik||X|| + ikaX3 — ioJmt) (7) 

for three different branches I = 1,2,3, where X|| and k|| are two-dimensional vectors with components (2:1,0:2,0) and 
(fci,fc2,0) = fc(cos#, sin#, 0) respectively (the angle 9 is counted from [100] direction), and = s^^Jk 2 + k 2 with 
the bulk sound velocity S(j) = S(/)(6*, </>), where cos</> = 

The solution for the phase velocity can be determined by substitution the plane wave of Eq. (Q) into the equations 
of motion ([!]). It gives the homogeneous set of linear equations. Setting the determinant of the coefficients equal to 
zero, produces a cubic equation in v 2 . Three roots of this equation are sq uares o f the velocities for three bulk phonons. 

For the propagation along (001) plane one of the velocities t\ = yj C44 / p is independent of 9 and represents a 
transverse mode. Two others depend on angle of propagations in the plane. They are neither longitudinal nor 
transverse, but we denote the upper branch by the letter "Z" and the lower one by H2" ■ 

The interface between two half-infinite media introduces an inhomogeneity in 2:3 direction. Therefore, we could 
expect that the plane waves also become inhomogeneous in this direction. The frequency and wave vector k|| are the 
same in both media, but the component fc 3 of wave vector may be complex and different in upper and lower media. 
Moreover, since boundary conditions Eqs.(|^^) comprise of 6 equations, the simplest solution for arbitrary direction 
of propagation should consist of linear superposition of three terms described by Eq. (Q) for each medium with their 
own different complex components k 3 . 

The further analysis is facilitated by performing a rotation of the coordinate frame in such a way that the direction 
of propagation of the acoustic wave in the plane of interface is along the x\ axis, i.e. k|| = (fc,0, 0). Let 



cos 9 sin 9 
— sin 9 cos 9 
1 



(8) 



be the transformation matrix which produces this rotation. Then, the transformation law for the elements of the 
elastic modulus tensor under this rotation is 

Kjki = ^ Tii'Tjj'Tkk'Tu'Xi'j'k'i'. (9) 

i'j'k'V 

At the first stage of the analysis we determine the possible values of k% for given magnitude of fc. To this end, we 
define fc 3 = ik[3, to — kv and we suggest the solution for the interface waves in the rotated system in the following 
form 

u.i = A i e~ kf3x ' i exp [ik(xx - vt)} , for x 3 > 0, 

u\ = A 1 ^'^ exp [ik{xi - vt)] , for a; 3 < 0. (10) 

Conceptually, the cc 3 -dependence is the part of the "amplitude" (See Ref.E3) and the wave-like properties are described 
by a common propagation part exp [ik{x\ — vt)]. Thus, the propagation vector is always assumed to be parallel to 
the interface even though the exponent (3 may be complex. 

If Rc/3, (3' > 0, then such a form describes a wave that propagates in x\ direction, whose amplitudes decays 
exponentially with increasing distance into the medium from interface. The waves with (i) Re/3 < and Im/3 < 0, (ii) 
Re/3' < and Im/3' < 0, (iii) Re/3, /3' < and Im/3, /3' < are the leaky waves which radiate the energy outward the 
interface. 

Substituting Eq. (|^Jl^) into the equations of motion Eq. (|l|) yields the set of homogeneous equations for each media 

L l0 {v,fi)A ] =0. 

^>,/3')'4=0, (11) 

where matrix L (or L') has the form 
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The sign plus (minus) corresponds to the upper (lower) medium, and we omitted primes for the lower medium, 
each medium Eqs. (O) have nontrivial solutions if the corresponding determinant of the coefficients vanishes: 

det(L) = 0, 



In 



(13) 

It gives the secular equation on unknown values of (3 with a phase velocity v as a parameter. The explicit form of Eq. 
( |l3| ) is given in Appendix |A| Due to the fact that we are seeking the solution for wave propagation in crystal plane 
of mirror symmetry this equation is bicubic in [3 and the roots have inversion symmetry with respect to the origin of 



One can also show that if /?/ = finj + iflij at j = 1, . . . , 6 are the roots of Eq. ( |l3| ) with complex 
ivi, then the roots f3j = /3r,- — iffij are the roots of the same equation with v = vr — ivj. Here the 



the complex plane 
velocity v = vr + 

superscripts "R" and "I" denote the real and imaginary parts respectively. 
The amplitudes A a (or A' a ) for any 0j {(3'j) are related by 



0) 



.4 



0') 



.4 



U) 



a 



a 



= K,, j = i, 



,6, 



(14) 



where the Kj are constants and Ca (v, (a = 1, 2, 3) are the cofactors of the elements in the first row of the matrix 
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The next step of our analysis is a construction of the general solution, which satisfies boundary conditions Eq. 
(II). To this end, we form a linear combinations from three terms ( J10[ ) with undefined constants Kj and -ftTj for each 
medium: 



3 r (j) 



i a = —jjrKj exp (ik(x\ — vt + i/3j(v)xs)) , for x% > 0, 
j=i c \ 
3 C 'U) 

i' a = 22 "(j) Kj ex P (ik{xi — vt — i/3j(v)x3)) , for x 3 > 0. 
j=i c \ 



(15) 



Substitution of this form for the displacement field into Eqs. (| ||) leads to a set of 6 (in general case) homogeneous 
linear equations for the Kj, Kj. The nontrivial solutions exist if the corresponding determinant vanishes: 



\D%\v)\=0, fc,l = l,. 



,6. 



(16) 



Eq. ( |l6| ) is the dispersion relation for the phase velocity v of the interface acoustic wave. In general it has to be solved 
numerically. The left hand side function D(v) = |-Djy | is some algebraic expression. Therefore, in general, the roots 
of Eq. ( |16|) are complex. Moreover, since D(v) comprises of six different decay constants /3j(v),/3j, which involve 
square roots from some expressions of v 2 , the function D(v) is multi-valued analytical function of complex variable 
v, defined on its associated Riemann sheets. The upper-script 7 enumerates these sheets. The number of Riemann 
sheets arc determined by different combinations of (3 branches. However, not all from 6!/3!3! = 20 combinations of j3j 
are possible for each medium in the superpositions of Eq. (Il5|) . Each of three decay constants must be taken from 
the different roots /3 2 of cubic equation det(L) = at fixed irT Therefore, the total number of possible combinations 
is 2 3 x 2 3 = 64. This is the number of Riemann sheets for our case. Since simultaneous change of all signs (3, (3' 
in Eq. ( |l6| ) does not change the form of determinant (see Appendix A), it is enough, in fact, to investigate 32 
independent Riemann sheets in order to find all possible roots of the dispersion relation. In the isotropic case and for 
the propagation along the directions of high symmetry the number of independent Riemann sheets reduces to 8. 
The sign convention for the sheets is determined by real part of (3. It is denoted as follows: 



(signRe(/3i), signRe(/3 2 ), signRe(^ 3 ), signRe(/3j ), signRe(/? 2 ), signRe^)) 



(17) 
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Let us assume that 7=1 corresponds to the case (++++++). If a solution exists on this sheet, then it is a true 
interface wave, which is also called Stoneley wave. All other 7 correspond either leaky waves or non-physical solutions. 
Some of them may also correspond to bulk phonons (see section III). In the next section we formulate selection rules 
for physical solutions. The right direction of the energy flux is the main principle for the selection. The explicit form 
of Eq. (16) and the way of enumeration of its Riemann sheets are given in Appendix 



III. SELECTION RULES FOR VELOCITY AND PHYSICAL MEANING OF LEAKY WAVES 

The question of total number of possible values of velocity has been investigated in earlier 70s theoreticallylll for 
the isotropic solid - liquid interface, and numerically for the case of two isotropic solidso. In the case of the liquid - 
solid interface there are eight Riemann sheets. It is shownt£l that the roots on all these sheets are the roots of an eight 
order polynomial in v 2 with real coefficients, and so there are eight complex roots which are either real or come in. 
complex conjugate pairs. Numerical investigation of the Stoneley equation ( A20| ) for isotropic solid - solid interfaced 



has shown that there are sixteen independent roots on its sixteen Riemann sheets. 

Thus, we can put forward a simple hypothesis: The number of possible values of v 2 is equal to the number of 
Riemann sheets. However, some of them may be degenerate so we are speaking about maximum number of different 
v 2 . Note that this hypothesis is true for isotropic surface wave either. It follows from it that the maximum number 
of modes in our interface is 64. To check this hypothesis we have calculated this number for one of direction which 
does not have any special symmetry The result is 64. 

Let us turn now to the problem of roots classification for the case 7 > 1. First of all we discuss real roots for v 
of the dispersion equation Eq. (|l6|). If all (3 are pure imaginary, it corresponds to refraction of bulk phonons and 
has nothing in common with leaky waves. If at least one of (3 at any side has a negative real part, such solution 
should be considered as non-physical. If it happens that some (3 are imaginary but some are complex with positive 
real part, then it relates to the problem of total internal reflection of bulk phonons. In numerical analysis we discard 
such solution, since they have a different nature. 

Now we come to complex v. Let us assume that we have real positive frequency 10 > and complex root 

v — v R =F ivi (18) 

of the Eq. (|l6|). As follows from Appendix A all complex roots form such pairs. Then, the wave vectors of propagation 
along xi-axis for these solutions will be complex and equal to 

k = lu/v = lu/(vr =F ivi) — hp ± ikj. (19) 

If kn > 0, it corresponds to the running wave propagating from the left to the right with exponentially decreasing or 
increasing amplitudes. Since both of these solutions always meet in pairs we will consider only the wave attenuated 
from the left to the right. Then we should chose only the root kR + ikj or vr — ivj, where vr, vi, kn 7 ki > 0. 
Since 7 > 1, one or several (3 have negative real part, i. e. 

13 = -p R - ipi, (20) 

where [3r > and sign /3j is not determined yet. 
It is useful to introduce the following notations: 

VR = (v 2 R + vj)/v R , (21) 
Pr = Pr ~ PiVi/vr, (22) 
Pi = Pi + PrVi/v r , (23) 

After substitution k — lu/v the term has a form of inhomogeneous plane wave (see also the succeeded section) 

e —uj(t—x/vR—z0i/vn) e —uj/vR(xvi/vR—zffii) ^24) 



In fact, the sign of Pi is not arbitrary, but is dictated by the radiation conditionliZllla. Indeed, since we consider 
here lossless media the only reason for the amplitude attenuation along direction of propagation on the interface is 
radiation of energy away from the interface into the bulk media. The total wave vector q = u>/vr(1, 0, Pi) is no longer 
parallel to the boundary, but is inclined to it, which indicates the presence of continuous flow of energy from the 
boundary to the bulk. Note that direction of phase propagation q/q does not represent the direction of energy flow 
itself. The later is determined by the time averaged power flux 
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W a = -Re{±ap a u* ), a,/3 = l,2,3. 



(25) 



One can show that for small vi the sign of the component W3 which determines the flow perpendicular to the interface 
coincides with the sign of f3i. Therefore, in notations of Eq. ( p0| ) the sign f3i must be positive to guarantee the proper 
radiation condition. n 

Now we come to the physical meaning of the leaky waves. Discussing Rayleigh waves Landau and LifshitzEl prescribe 
to drop as non-physical the solution, which increases away from the surface. In the theory of leaky waves we take into 
consideration such a solution. The goal of this part is to give physical explanation of leaky waves (see also review by 
MaradudinB). 

On both side of the interface our solution consists of three inhomogeneous plane waves. At large values of x and 
z these waves can be considered as independent. So we concentrate on one of them choosing the wave at z > with 
negative real part of /3, i.e. the wave exponentially increasing with z. This wave has a form Eq. (|2^). 

e ~iu(t-x/v R -zi3 r /v R ^ e -uj/v n (xvr /v n -z/3 R ) (26) 

where (3r > and /?/ > 0. This is an inhomogeneous bulk plane wave with wave vector q = lo/vr(1, 0, /3j) propagating 
away from the interface. The lines of constant phase are determined by the equation t — x/cr — zfii /cr = Ci, and 
the lines of constant amplitude are given by the equation xvj /vr — z(3r — C-i (See Fig. 2) The later expression shows 
that such modes attenuate when they move along the surface z = 0. 




V R 

FIG. 2. Illustration for geometry of leaky interface wave term. Power flux W and wave vector q are shown. The full lines 
with different thickness are the lines of a constant amplitude. The thickness of the lines indicates schematically the absolute 
value of the amplitudes. Dashed lines are the lines of a constant phase. For visual clarity the angle of constant amplitude lines 
has been exaggerated. 



The central point of our understanding of the leaky waves is that one can not consider this waves in whole region 
of x, since amplitude becomes infinite when x — > — oo. This happens because we have chosen the solution which 
propagates from the left to the right along cc-axis. Thus, we should propose that the wave creates at some line, say 
x = 0, in the plane of the interface. Our equations of motion do not include any dissipation, therefore the attenuation 
in the plane of interface may be only due to radiation into the media. There is an important theoremEil stating that 
for inhomogeneous plane waves the energy flux is parallel to the plane of constant amplitude. The cross-sections of 
these planes with the plane zx are shown by full lines of different thickness. The upper line is thicker because the 
amplitude of the wave at the line x = is the largest and it decrease with increasing x. That is why the amplitude 
at any point x increases with z at z < xvij ' (3vr. It follows from the above theorem that the energy flux can not cross 
the planes of constant amplitude. It also can not cross the plane of maximal amplitude. This means that all the wave 
is within the wedge formed by the plane of interface and the plane of maximal amplitude z = xvj / PrVr, at least in a 
sense that the whole energy of the wave is within this wedge. The amplitude of the wave is finite everywhere in this 
region. 

One gets a severe contradiction considering a stationary problem. Indeed, the equation (^6|) gives nonzero result 
outside the wedge as well. Moreover amplitude diverges when z tends to infinity. This is an artifact of the stationary 
consideration. The origin of the divergence is the infinite amplitude of the wave at point x = — oo. The increase of 
the amplitude at large z is an artifact originated from the flux coming from large negative x. 
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It is important to mention that the leaky wave does not differ substantially from the true interface wave only if 
kji 3> kj in Eq. (18) or vr vj. Since (3r is not small, it means that the angle between the planes, forming the wage 
in Fig. 2, should be small. If this condition is not satisfied, the interface (or surface ) wave can not be considered as 
a wave since the wavelength is larger than the attenuation length. 

The problem does not contain any small parameter which could make this condition fulfilled. The numerics show, 
however, that the majority of modes have small attenuation. The roots of this phenomenon are not clear for us. 

Thus, based on discussion in this section we use the following selection rules for values v — vr — ivi, which are 
solutions of Eq. ( |l6| ) . 

1. v R > 0, vi > 0. 

2. If vi = 0, than Re/3, Re/3' > 0. 

3. If vi > and Re/3 < 0, than Im/3 < 0. 

4. If vi > and Re/3' < 0, than Im/3' < 0. 



IV. NUMERICAL RESULTS 



To calculate velocities as a function of angle we use Eq. (Al) for (3 2 and Eq. (A6). At the first step we divide 
complex plane v in interval 2 km/s< vr < 7 km/s, < Vi < 0.6 km/s into 400 squares. For t he v ertex of each square 
we find six values of (3 for upper medium and six values of (3' for the lower medium using Eq. ( Al). For each value of v 
we find 32 different combinations of j3 (each of 6 values) and substitute them into Eq. ( jA6| ) . For each Riemann sheet 
7 we find all minima of abs(|-D^(«)|) with respect to vr and Uj using a standard program. For those minima which 
are close to zero, we do iterative search of roots. The parameters of the bulk lattices have been taken from the Table I. 

ities(g/cm 3 ), elastic constants (10 10 N/m 2 ), and sound velocities (km/s) in the directions of high symmetry 



TABLE I. Dm 
for bulk crystals!! 



Crystal 


P 


Cn 


C12 


C44 


'[100] 


ti 


'[110] 


ta[no] 


Alo.3Gao.7As 


4.794 


12.24 


5.65 


5.90 


5.05 


3.51 


5.56 


2.62 


GaAs 


5.307 


12.26 


5.71 


6.0 


4.81 


3.36 


5.31 


2.48 



To check-thc method some results have been obtained using completely different Surface Green Function Matching 
methodE3o, which is discussed in Appendix B. We have not found any differences between the results of two methods. 

For an additional check of the problem we have calculated true surface waves for both materials. The dispersion 
relation for surface waves on (001) cut of cubic crystals can be obtained from the determinant of truncated interface 
matrix Eq. ( |A5| ) . For the upper (lower) medium we should take lower left (right) 3x3 part of the matrix. We hawft 
gotten 2.873 km/s (2.737 km/s) for GaAs (Alo.3Gao.7As). These results may be compared with the results by Farnellc£l 
Our results are slightly different because of the difference in the parameters of bulk materials. Taking parameters 
used by Farnell we have obtained his velocities with a very high accuracy. 

The results for leaky interface waves are shown in Fig. |^ and Fig. ^. Plotting this figure we have taken into account 
all selection rules formulated at the end of section III . The discontinuities appear because at some points these rules 
are not fulfilled and the corresponding modes become non-physical. All the modes have different anisotropy, different 
attenuation and different angle intervals of their existence. In Fig ^3] we draw all the modes (Note that velocity scales 
in 3(a) and 3(b) are different). One can see that there are abundant number of leaky interface modes in the velocity 
range between 3km/s and 4km/s. Thorough analysis shows however that a majority of those modes exists in a very 
restricted range of the angles. Probably, it will be difficult to detect them experimentally. In Figure f| we draw 
selected modes which exist in the large interval of the angles only. We left one mode (A) with a small angle range in 
the figure as a typical example of discarded modes. For further clarification of the picture we also discard all physical 
modes with a strong anisotropy of their real and/or imaginary parts of the velocity. One of such modes with strong 
real part anisotropy (B) is left in the figure as an example. 



7 




FIG. 3. The real parts of complex velocities for all leaky interface waves in the range 3.5 km/s < vr < 5.7 km/s (a) and 2.4 
km/s < Vr < 3.5 km/s (b). The solid lines present velocities of bulk acoustic waves for both media. 
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FIG. 4. The real parts of complex velocities for the selected modes in the range 2.4 km/s < vr < 5.7 km/s. See explanation 
in the text. The numerical parameters for these modes are given in the Table II. The value of attenuation is determined by the 
imaginary part of complex velocity. Different linestyles and thickness correspond to different magnitude of average imaginary 
part W (km/s) of velocity. The solid lines present velocities of bulk acoustic waves and the dotted lines are true surface acoustic 
waves for both media (Upper curves correspond to Alo.3Gao.7As and lower curves correspond to GaAls). 



All numerical parameters for the modes presented in Fig. 
introduce the following parameters. 



are given in the table II. For modes classification we 







1. Averaged over angles real vr and imaginary vj part of the velocity. 

2. Angle range parameter S an — (tp m a.x ~ fmin) 4/7r. It equals to unity when a mode exists on whole range of angles 
and less than one otherwise. 

3. Anisotropy parameters for real and imaginary parts of the velocity 

0»H = (fKmax - WflminV^fl, Vl = {Vj max ~ W/min)/W- (27) 

Modes with smaller anisotropy have smaller a. 

Among the modes there are two ( C and D m, the figure) which remind those in the surface acoustic problem. Namely, 
for Al^Gai-^As materials there are alwaysc3 true SAWs (they are shown by dotted lines in the figure), which change 
very a little with angle until they meet the lowest bulk transverse velocity curves (they are shown by solid line) . After 
that point the true SAWs repeat behavior of their bulk velocity curves up to the end - 45 degrees angle on the picture. 
Meanwhile, in the region between two transverse bulk velocity curves the leaky surface acoustic waves appear with 
approximately the same real part of velocity. 

The situation for interface waves is different. As we mentioned, there are no true interface waves. However, there 
are leaky interface waves with a very small attenuation at small angles before "colliding" with the bulk transverse 
velocity curve. At larger angles the leaky modes acquire larger imaginary part of the velocity and become stronger 
attenuated (See Table II). 

V. CONCLUSION 



We have derived a dispersion equation for interface waves at the interface of two cubic crystals in the plane, 
perpendicular to [001]. Analyzing different solutions for the interface waves we come to conclusion that the total 
amount of interface modes in each direction is equal to the number of the Riemann sheets. In our case this is 64. We 
have successfully checked this hypothesis by calculating the number of modes in one direction of the interface plane, 
which does not have any special symmetry. 

The computations have been made for the interface Alo.3Gao.7As/GaAs. We have shown that in this case all 
interface modes, which have physical meaning are leaky, but majority of them have small attenuation in the direction 
of propagation. We show that for understanding of the physical meaning of leaky waves one should consider not a 
stationary problem , but the problem starting with creation of the wave at some line in the interface plane. 

After that we are able to formulate how to separate this 64 modes into physical and non-physical modes. This 
separation mainly based upon some theorems on the energy flux and upon an assumption that if a mode deviates 
from the interface in some medium, the energy flux should go in the same medium. 

Using the elastic moduli of the bulk lattices we have performed numerical calculations of the velocities of the 
interface waves as a function of an angle in the plane of interface in a wide range of velocities. The results are shown 
in Fig. [|. One can see two close groups of modes within the intervals 3—3.5 km/s and 4.2—4.5 km/s respectively. 
These groups may be responsible for two periods of oscillations which have been observed in the experiment with the 
two-dimensional electron gas in magnetic fieldlij, mentioned in the introduction. Note that the velocities of the leaky 
interface waves may be sensitive to the difference of the bulk media parameters. This difference is not known good 
enough. This fact may be responsible for possible deviation of our calculations from the experimental data. 
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TABLE II. Parameters for the selected modes. The units for average velocities are km/s. 
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APPENDIX A: EXPLICIT FORM OF THE DISPERSION RELATION IN GENERAL, SYMMETRY 

DIRECTIONS, AND ISOTROPIC CASES 



From the determinant (13) we obtain the following equation on unknown variable (3 with a phase velocity v as a 
parameter: 



/3 6 -/3 4 



a +b + c — 



(A 2 - 1 - df 



A 2 



+ /? 2 



ab + bc + ca — b 



(A 2 - 1 - d) 2 
A 2 



where 



= A 2 



1 - 

2 



b = 



1 



1 

t = - sin AB. 
4 

^'Weight factors" in Eq. (|l5| ) have the following form: 

1, 



v 

T 2 



(J) 



L21 

L22 

_L 3 i - 



g/4sin4fl 



-C 4 4/3|+C 4 4+-D/2s 



Ci 1 — C4 



i — pv* 
-D 



+ c(t 2 - ab) = 0, 



(Al) 



for a = 1, 
for a = 2, 

for a = 3, 



(A2) 



(A3) 



where the upper sign is for the upper medium, and for a convenience sake we omitted primes in formulae for the 
lower medium. Then, for each term in the sum ([l5]) the condition (Jfy on the stress continuity in rotated frame can 
be written in the vector form: 
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(A4) 



One can can see that the boundary conditions decouple for the sagittal plane {xix^} and the perpendicular direction 
X2- With the help of Eq. (A.4) and minor simplifications we obtain six homogeneous linear equations for the boundary 
conditions with the matrix L)(v) equals to 



C 4 4/3i(1-P^) 
C12 + Cufipl 1 



pf 
fopf 

C44/3 2 (1 - P i 2) ) 



C44/3 3 (l-pi 3) ) 

C 12 + C u plpf 5 



-1 
'(1) 

-J>2 

/3;^ (1) 

C 4 4/3ip 2 (1) 
-C12 - 



-1 

- P2 (2) 

/? 2 % (2) 

12 — ^llP2 P3 



-p' 2 (3) 



4' 



-a 



(A5) 



The condition for a nontrivial solution is that its determinant must vanish: 



D^(v;P u ...,P' 3 )^\D\]>\=0 



n - ln(7)i 



(A6) 



Since each (3 is doubled value function of complex variable v, this determinant has 64 Riemann sheets, which we 
denote by upper-script 7. 

The sign convention for independent sheets is determined by real part of (3. It is given as follows: 



sign (Re(/3i),Re(/3 2 ), Re(/3 3 ), Re(^), Re(#),Re(&), ) . 



(A7) 



That is, Riemann sheet (+ + 



-) corresponds to the case Re(/3j) > and Re (/?'■) < 0, where j = 1, 2, 3. 



The important property of the determinant is that its solutions are invariant against the simultaneous change of 
sign of all six decay constants. Indee d, t he functions P2, P3 are defined in such a way that they depend on j3 2 only. 
Therefore, from the form of matrix ( |A5| ) it follows that simultaneous change of sign (3 leaves the secular equation 
unaltered. Thus, in order to find all the roots of Eq. ( [A6] ) it is en ough to investigate only 32 independent Riemann 
sheets. To enumerate these sheets we, firstly, solve equation (Al) and sort obtained (3(v) for each medium in the 
following order: 

|Re(/?x)| < |Re(/? 2 )| < |Re(&)|. 

Secondly, we choose notation that the sign of the real part f3\ is always positive. Then the upper Riemann sheet 
7 = 1 corresponds to the case (++++++) and the subsequent numbers for lower sheets are given in the table III. 
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Another feature of Eq. (A6) is the following. If v = vr + ivj is the solution of the dispersion relation D^'(v) = 



Dfi^+iDf' ~ 0, then v = v^ — ivi will be also a solution of Eq. (|A6|). Here we explicitly separated real and imaginary 



,(7) _ 
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part of complex function. It could be understood by the following consideration. Let 0i = 0r 1 + i(3i [ (I = 1, . . . , 6), 
obtained from Eqs. ( |l3|) for v — vr + ivj, then the solutions of Eqs. ( [l3l) for v = vr — ivj are 0\ = fin, — i0i r 
Moreover, the simultaneous change v — > Vr — ivj, 0\ — > 0r 1 — i/3j l in Eq. (A6) brings up the change of sign at the 
imaginary part of the determinant only: 



D^\v R - ivi) = D 



(7) 



iD 



(7) 



Thus v = Vr — iv i is also the solution. 

If we suggest that frequency of the interface acoustic wave is real, then k — lo/v will be complex. We will choose 
solutions with v = Vr — ivi which correspond to the attenuated waves (for u> > 0) along direction of its propagation 
on the interface. 

The dispersion relation simplifies considerably for the propagation on the interface (001) along the directions of 
high symmetry ([100] and [110]) and in the isotropic case. For all this cases, matrix elements L12 = L21 = [because 
of sin 40 = or D = 0, see Eq. (|l2|)]. Therefore, the systems ( pi] ) for each media breaks up into the pair: 
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L13 £33 




^3 



= 



(A8) 



and 



LonA-2 



0. 



(A9) 



The only nontrivial solutioii|-ef the latter corresponds to a bulk transverse acoustic wave propagating parallel to the 
interface of the elastic mediaEJ, and therefore it is discarded. Thus, in these cases the interface waves polarized in the 
sagittal plane and do not have ui component. 

For the interface waves polarized in the sagittal plane the solvability condition for (AS) is the biquadratic equation 



A 2 /? 4 - (3 2 {ll + A 4 7 i - (A 2 - 1 - df) + A 2 7 2 7 f = 0, 



(A10) 



where we introduced notation 
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In the isotropic case (d = 0) the Eq. ( AlC| ) reduces to the 
with obvious solutions 



for 9 = 0, 
for 6 = tt/4. 
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(A15) 



(A16) 



Pi = ±^/l-v 2 /t 2 



2 = ±^l-v 2 /(\t) 2 



(A17) 



For all these cases the general solution for the interface wave are linear combinations of two partial waves for both 
media: 
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After substitution solution (A18) into the boundary conditions (@,§|) we obtain the set of four homogeneous linear 
equations with solvability condition: 



fhi& hvf P[pT ] &p? ] 

cwMi-p^) cMi-pf) c'Mi-p'P) c'Mi-p'P) 



(A19) 



For isotropic case p^ 1 — — 1 / /3f and p { 3 ' = — 1. Then the resulting dispersion relation reduced 2 ! to 

v* [(p - p'f - (p(3' 2 + p'0 2 )(p0[ + p'0i)] + 4iV [p(3[(3' 2 - p' (3,(32 -p + P '}+ 4F 2 (1 - /3 1/ 3 2 )(1 - (3[f3' 2 ) = 0, (A20) 
where 



,(2) 



F = pt — p'f 



1 ,12 



(A21) 



There are sixteen independent roots on eight different Riemann sheets for the given equation There is always non- 
attenuated solution for t — t' , A = A' and p ^ p' . 



APPENDIX B: GFSM METHOD 



In addition to our method described above it is also possible to obtain the dispersion relation for the interface wave 
using Surface Green Function Matching (SGFM) analysis (see for details Refs.EjO). We used this technique as a 
complementary independent method in our numerical calculations to compare the results. The main advantage of the 
method is that it gives the determinant of the matrix with dimensions 3x3, rather than 6x6. However, the calculation 
of the those matrix elements is more cumbersome. Another important development of this method is construction of 
the Green function of whole system, which also enables to calculate interface contribution to the important physical 
properties of the system such as the density of phonon modes, specific heat, and the atom mean square displacements. 

We consider the (001) interface of cubic crystals and arbitrary propagation directions. The idea of the method 
is to construct the Green function G s of the composite two-media system, such that it incorporates the boundary 
conditions at the interface. In fact, to obtain the dispersion relation it suffices to know the surface projection g s of 
G s on the interface plane. Since the boundary conditions (|B|,H) involve first derivatives, the surface projection of the 
bulk Green functions Gij and their normal derivatives will be involved. These are defined by 
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Gij(q)iqexp(=piqTi) dq. 



In our case G(q) = ft] , where matrix L is given by Eq. ( |l2| ) with substitution i(3 — > q. The calculation is tedious 
but straightforward 2 ^. 
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Here 6* is the angle formed by the direction of propagation and the Xi axis, r = — s "^ 4g and notations for A,d,7i,t 
are given by Eqs. (All- A14). The pi, /3 2 and /3a are the solutions of Eq. (p"3|). It is provedES that g s is given by 

^-^aWg-'-A'Hg'" 1 , (B2) 

where the inverse is defined in th e tw o-dimensional space of the interface. The matrix A is a suitably defined linear 
operator. Its appearance in Eq. (B2) reflects the boundary conditions. 
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The secular equation, which gives the interface mode dispersion relation is shown to be 

det [g 8 ]- 1 = 



(B3) 



(B4) 



For high-symmetry directions the further simplifications are possible. In this case we separate the transverse bulk 
mode again and the resulting matrix [g s ] _1 for interface waves is two-dimensional with the following matrix elements: 
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G 4 4A 2 /3 1 /3 2 (/3 1 2 -/3 2 ) 
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Pi and fa are the solutions of Eq. (A10) with 7 2 defined by Eq. (A15). The dispersion relation for the interface 
waves can be obtained from the zeros of det [6/ 17 s ] 

The former expression for the determinant can be reduced in a very easy way to the relationships for the isotropic 
case. When D — Cu — C12 — 2G 44 = the matrix elements are given by: 
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